Spectrum and eigen functions of the operator 

Hjjf(x) = f(U — l/x)/x 2 and strange attractor's density for the 

mapping x n+1 = 1/(U - x n ) 

Kozlov G.G. 

We solve exactly the spectral problem for the non-Hermitian operator Hjjf(x) = 
f(U — l/x)/x 2 . Despite the absence of orthogonality, the eigen functions of this 
operator allow one to construct in a simple way the expansion of an arbitrary function 
in series. Explicit formulas for the expansion coefficients are presented. This problem 
is shown to be connected with that of calculating the strange attractor's density for 
the map x n +i = 1/(U — x n ). The explicit formula for the strange attractor's density 
for this map is derived. All results are confirmed by direct computer simulations. 



The author of this note is not an expert in the theory of strange attractors. So, the 
results presented below may appear to be well-known to those who work in this field. If 
the material of this note will be found familiar, totally or to some extend, for some of the 
readers, the author will much appreciate sending him appropriate references which will be 
included in to the next version of this publication. 

Let us consider the linear operator Hjj, depending on the parameter U and acting on the 
arbitrary function f(x) (with real x) in accordance with the following definition 

Huf{t)s ^f(y-1.) (1) 

and consider the spectral problem for this operator i.e. find its eigen numbers \y and eigen 
functions &u(x): 

H u( j$j{x) = \u<Tu(x) (2) 

We solve this problem under the following conditions: (i) the eigen functions of (1) must 
be finite, (ii) the parameter U obeys the inequality \U\ < 2 which corresponds to the regime 
of strange attractor for the mapping x — > U — 1/x which is co inside with operator (1) up to 
the factor 1/x 2 providing the equality of the integrals / f(x)dx = J Hjjf(x)dx. 

Using (2) and (1) we write (we omit indexes U and A): 

a{x) = \L AU _ 1/x) = ^_L_f a{u _ m _ l/x) = (3) 
T' Uu-i/w-i/u-i/J *" - 1/u - 1/u - 1/x) = 

n 

, a „ n 

U x (u-i/ x){ u-i/u-iM... h u - 1 < U - l ' u - l l- l ' x > 

Define T as: 

Tx = U - 1/x (4) 
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Then one can rewrite (3) as: 



r n— 1 



a{x) = 



X' 



1 -2 



k=0 



X 



a(T n x) 



(5) 



Now note that the eigen function of (1) Hjj with the eigen number equal to unit A = 1 can 
be written in the form: 



<7$j=\x) = Cu(x) EE 



1 



2tt x 2 - Ux + 1 



(6) 



We normalize this function by the condition / Cjj(x)dx = 1. 
Now one can express the complex fractions entering (3) as: 



rTl-l 

]jT k x 

k=0 



1 -2 



C v {T n x) 



and consequently 



a(X 



Crr(x) 



(7) 



(8) 



a(T n x) X n Cu(T n x) 

This formula is valid for an arbitrary eigen function of (1) and for any integer n. Taking into 
account the fmiteness of the eigen functions of (1) we come to the conclusion that |A| = 1, 
Otherwise, the left-hand side of (8) either vanishes or diverges at n — ■> oo which contradicts 
to the finiteness of the eigen functions a(x). Thus: 

A = exp(«0) (9) 

with real 0. Equalising the module of both parts of (8) we obtain: 

\a(x)\=Cu(x) (10) 

and consequently: 

a(x) = Cjj{x) exp(i9(x)), (11) 
with the function 6{x) being real. Using (8) and (9) we obtain for 9(x) the following equation: 

expt[n<f)-6(T n x) +6(x)\ = 1 (12) 

whence 

n<f) - 6(T n x) + 6{x) = 27rMi n (a;) (13) 

where Mi n (x) is an arbitrary function taking the integer values. Suppose we find the 
function 9{x) obeying: 

6{T x) - 6(x) = 6{U - 1/x) - 6(x) = + 2ttM(x) 

with M(x) being an integer function of x. Then: 

6(T k+1 x) - 6(T k x) = <f) + 2ixM{T k x). 



(14) 



(15) 
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Now summing (15) over k from zero to n — 1 we obtain: 



n— 1 



6{T n x) -9{x) =n0 + 2vr^M(T fc x). (16) 



k=0 



Here T°a; = x. Consequently the function 6{x) defined from (14) satisfy (13) with 



n-l 



Mi n (x) = -J2 M(T k x) 

k=0 

Thus the problem reduced to solving the equation (14). We will see below that solution of 
this equation can be found only for some particular value of angle 0. 
To solve the equation (14) let us consider an auxiliary equation 

G(x) ? 1 ; 

Note that hnG(x) obey the equation similar to (14). Construct the solution of (17) in 
the form: 

v ' 1 + bx 

Substituting (18) in to (17) one can obtain the set of equations for a, b £. The solution 
of this set of equations gives: 



C _ U*-2±*\U\V4=U* (ig) 
- 1 + { (20) 



u 

» = -i^ (21) 

If \U\ < 2 (the strange attractor regime for mapping (4)) then |£| = 1. The sign of square 
root in (19) is not important. Then: 

£ = exp* 7 (22) 

with 

ry - aic UB^IEE.) + J when V2<\U\<2 , . 

7_arctg^ u2 _ 2 ) + ^ whenV 2 > 1^ ^) 

Here we take into account the continuity of £(U) function (19). Then it is easy to see that 
|G(x)| = 1 and consequently: 

G{x) = expzae (24) 



with 



ae xsm7 
tg 2 = z(l + cos 7 )-tf { ' 



and consequently: 

x sin 7 



ae = 2 arctg 



x(l + cos 7) — U 



+ 2irM(x) (26) 
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as: 

.c sin o 

9(x) = 2 arctg 



where M(x) is some integer function of x. Then one can see that 

ImlnG'(x) =ae(x) (27) 
satisfy the equation (14) with = 7- Thus the final solution of equation (14) can be written 

x(l + cos0) - u\ ^ 

f \UWA-U 2 \ Jo when V2 < \U\ < 2 

= arCtg l IP -2 ) + 1 n when ^2 > (29) 

and 

6(U-l/x)-6(x) = <j> + 2nM(x) (30) 

Now one can see that the function 9 n (x) = n 9(x) where n is integer, satisfy the equation 
(14) with — > <f) n — nip. 

Thus one can write the set of eigen functions <Ju(x) and eigen numbers A n of (1) as: 

(Jij(x) = Cu(x) expm9(x), \ n = expzn(p (31) 

where n is integer, and £jj(x), 0{x) and are defined by (6), (28) and (29). 

I. PROPERTIES OF THE EIGEN FUNCTIONS AND THE RELEVANT 

EXPANSION 

Functions (31) should possess the following property: each of these functions (except for 
Cu(x) with A = 1) has zero integral. To show this let us integrate over x the equation for 
an arbitrary function 

/ % o{ V ~ l l x ) = X J °( x ) dx - ( 32 ) 
Replacing 1/x — > x in the left-hand side we obtain 

(1 - A) J a(x)dx = (33) 

Thus if A ^ 1, then 

J (T x ^{x)dx = (34) 
The next property of functions (31) which we use below is: 

<t%(x) exptm9u(x) = al +m (x) (35) 

These properties allows one to expand an arbitrary function f(x) in series using the set 
of functions (31) as follows. Suppose that function f(x) can be expanded in series 

+oo +oo 

f(x)= E C n aZ(x)= £ C n Cu(x)e^n0u(x)]. (36) 
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After multiplying both sides of (36) by exp —im6u{x) and integrating, one can obtain the 
following expression for the coefficients C m : 



J f(x) exp [ — im6u(x)]dx 
J Cjj(x)dx 



f(x) exp [ — im0u(x)]dx 



(37) 



The function 9jj(x) is discontinuous at 



x 



U 



1 + COS ( 



It can be converted into continuous one by adding proper function in the form 27rM(x) (with 
M(x) being integer function). The function 6u{x) defined in this way will be as follows: 



Ou{x) = 2 arctg 



xsmcf) 



x(l + COS0) — U 



+ 



when x < Xq 

—2n\U\/U when x > x 



(38) 



The possibility of expansion (36), (37) was verified by direct computer calculations for 
various functions f(x). 

By simple algebra functions (31) (Ju(x) can be rewritten as: 



a^(x) = C v {x) 



Cu(x) 



exp[i9u(x)] 



x(\U\ +1VA-U 2 ) -2U/\U\ V 
x(\U\ - iVI^TP) - 2U/\U\) 

VI^TP 1 



(39) 



2tt 



Ux + 1 



x{\U\ + iV4 - U 2 ) - 2U/\U\ 
x(\U\ - tV^TP) - 2U/\U\ 
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Omitting some unimportant coefficients one can finally present the eigen functions of (1) 
in the form: 

UK ' 2mU\x-Ru x-RfjJKx-RuJ V ; 

Ru = ^7 (42) 

RuRu = 1 (43) 
An expansion of an arbitrary function by functions (40) has the form: 

+oo 

f{x) = £ C n aZ(x) (44) 

n=— oo 

c »=JH^Y dx (45) 

The example of above expansion for various number of "harmonics" is presented on Fig.l. 



II. LINEAR-FRACTION MAPPING, AN EXPLICIT FORMULA FOR AN 
ARBITRARY NUMBER OF STEPS, STRANGE ATTRACTOR, CYCLES 



Consider an arbitrary linear-fractional mapping of the form: 

1 + ax n 

Xn + l = ~r~, • 

+ cx n 

One can reduce (46) to the form: 



(46) 



x "» = whr n (47) 

by replacing x n — > k\X n + k 2 with appropriate constants k\ k 2 - For this reason without loss 
of generality one can study mapping (47). Introducing new variable y n as 



x n = j^fL (48) 



and back 

Ox— yy 



1 — 9x 

one can reduce the corresponding mapping for variable y n to the simple form: 

Vn+i = aey n (50) 
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with the following values of parameters entering (48), (49) and (50): 



U - VU 2 - 4 



(51) 



or 



9 = 

33 

e = 
e = 



C/ + Vt/ 2 - 4 
t/ - Vt/ 2 - 4 

u + VTP^t 



U + Vf/ 2 - 4 



(52) 



33 = 



C/ - VC/ 2 - 4 
U + Vf/ 2 - 4 



[/ - Vf/ 2 - 4 

Thus if we perform mapping (47) n times using x as a starting point we obtain for x n 
the following expression: 



33 n y 



where y 



1 — 6xq 



(53) 



^ + eae"?/o JU ^0-1 

In the case of strange attractor regime (i.e. when U 2 — 4 < ) 33 = exp?0 with <p = 
2arctg[V4 — U 2 /U]. If = 2n/n, with n integer, then a; m+n = x m and we obtain the simple 
cyclic regime. For given n one can find the value of U for which the simple cyclic regime is 
appeared: 

A 

(54) 



U 2 = 



tg 2 (7r/n) + 1 



The general cyclic regime corresponds to the condition: ri(p = 2rrm where m and n are 
integer. 



III. DENSITY OF STRANGE ATTRACTOR FOR MAPPING (47) 

The values x n obtained in the process of successive mapping (47) are aperiodically walking 
along the axis of numbers (when \U\ < 2, strange attractor regime). Let us define the relative 
density of values x n after N steps as 

1 N 

JV n=0 

We are interested in the ultimate density function defined as 

1 N 



n— >oo M 

JV n=0 
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For this function one can obtain the functional equation as follows. Using the fact that 
all quantities x n are successively derived from each other by mapping (47) we can write: 



I N-\ f 1 \ I 

M = s E <* - —) + ^ - *>)■ <») 

Taking into account that 

'('-rb)-^— i) <67) 

we obtain: 

1^1/1 \ 1 1 

^) = J™ - r ^) + ^ - *o) - - 1/* - xn). (58) 



n=0 

Now one can see that at N — > oo: 



p(z) = \p(U-l/z). (59) 

,2- 



We see that the density we are interested in is the eigen function of the operator Hu, with 
eigen value equal to unit and therefore has the form of lorentzian (6) : 

J A - U 2 1 

Fig. 2 shows the density of strange attractor obtained by computer simulation for various 
values of parameter U. The smooth curves were calculated by (60). 

All the aforesaid is valid for the case of the strange attractor regime \U\ < 2, when 
the density is not trivial. When \U\ > 2, the density has a delta- singularity at z — 1/0 for 
|ae| < 1 and z — l/£ for |ae| > 1. Equation (59) can be generalized for the case of an arbitrary 
mapping of the form x n+ \ = f(x n ). Denoting the inverse mapping by F (i.e.x n = F(x n+ i), 
we obtain that the density of attractor of mapping / satisfies the equation: 



p(z) = F\z) P [f{z)) 



(61) 
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Captions 

Fig.l The example of expansion in to series using the set of eigen functions of non- 
hermitian operator (1) 

Fig. 2 The density of strange attractor obtained by computer simulation for various values 
of parameter U. The smooth curves were calculated by (60). 



